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SUMMARY 

In low cycle fatigue, cracks are initiated and propagated under 
general-yielding cyclic-loading. For general-yielding cyclic-loading, 
Dowling and Begley have shown that fatigue crack growth rate corre- 
lates well with the measured AJ. The correlation of da/dN with AJ 
has also been studied by a number of other investigators. However, 
none of these studies has correlated da/dN with AJ calculated speci- 
fically for the test specimens. 

Solomon measured fatigue crack growth in specimens in general- 
yielding cyclic-loading. The crack tip fields for Solomon's speci- 
mens are calculated using the finite element method and the J-values 
of Solomon's tests are evaluated. The measured crack growth rate 
in Solomon's specimens correlates very well with the calculated AJ. 

INTRODUCTION 

In strain-controlled fatigue or low cycle fatigue, cracks are 
initiated and propagated under general-yielding cyclic-loading. There- 
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fore, linear elastic fracture mechanics can no longer be used to analyze 
fatigue crack growth rate. 

Dowling and Begley [1] and Dowling [2] correlated fatigue crack 
growth rate with AJ as shown in Figure 1. The value of J was evaluated 
as the rate of change of the deformation work density with respect to 
crack extension for through-cracks in plates. The data points in the 
figure were measured with specimens under general-yielding cyclic-loading, 
and the scatter band was obtained with the linear elastic fracture mech- 
anics specimens. They agree very well with each other. Dowling [2] has 
developed an equation for calculating AJ for surface cracks by combining 
the finite element method calculations of Shih and Hutchinson [3] for 
center-cracked and edge-cracked panels. Haddad and Mukherjee [4] and 
Tanaka, Hoshide, and Nakata [5] followed the same procedure to evaluate 
J and correlated J with da/dN. Kaisand and Mowbray [6] correlated fatigue 
crack growth in general yielding with AJ. They divided AJ into two parts: 
elastic and plastic. 


AJ = AJ + AJ (1) 

e P 

Using Shih and Hutchinson's calculation and following a procedure similar 
to Dowling's, AJp for a surface crack is approximated by 

AJ =1.96 /1/n' AW a . (2) 

P P 


where n' = the cyclic strain hardening exponent and AW^ = the applied 
plastic deformation work density. 

Tomkins [7] and Tomkins , -Sumner , and Wareing [8] correlated crack 
growth rate with crack tip opening displacement, CTOD, and J. Here, J 


consists of two parts 
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2 airoe a 

j = + L_ (3) 

E 1+n 

where n is the monotonic strain hardening exponent. However, in all of 
the above investigations (1, 2, 4, 5, 6, 7 and 8), no attempt was made 
to correlate da/dN data with J calculated for the specific test specimens. 

Solomon [9] measured fatigue crack growth in specimens in general- 
yielding cyclic-loading . In this paper, the crack tip field of Solomon's 
specimens will be calculated with FEM. The J values will be evaluated, 
and the measured crack growth rate will be correlated with the calculated 
J values. 


ANALYSIS OF SOLOMON’S EXPERIMENT 

Solomon measured fatigue crack growth in 1018 steel under general- 

yielding cyclic-loading. The chemical composition of the steel is given 

in Table 1. The test section of the cylindrical specimen was reduced 

by two semi-curcular notches (Figure 2a). The gross cross-sectional area 
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of the test section of the specimens was 1.24 x 12.6 mm . A sharp notch 
was used to initiate the fatigue crack. Both ends of a specimen were 
rigidly attached to the loading frame of the test machine. The fatigue 
crack growth rates were measured at six different applied cyclic plastic 
strain ranges. The applied strain range was controlled by the extenso- 
meter located 7.62 mm from the edge of the specimen as shown in Figure 2a 
The total applied strain range consists of two parts 

Ac = Ac + Ac (4) 

e p 


Ac 


e 


and 


Ao 

E 


(5) 
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Solomon [9] plotted his crack growth data in terms of pseudo stress- 
intensity factor defined as 

A(PK) = E(Ae) /a (6) 

Subsequently, the data were analyzed by Haigh and Skelton [10] in 
terms of a strain-intensity factor defined as 

AK = (^ Ac + Ac ) Aa (7) 

c. z e p 

For comparison with (7), (6) can be written in the form 

^ - (P ^ = (Ac + Ac ) AI (8) 

E e p ' 

Solomon correlated his data with the total strain range Ac. In Haigh 
and Skelton's analysis, the elastic strain amplitude was used instead of 
the elastic strain range because of the consideration of crack closure 
when the applied stress became compression. The data and the correlations 
proposed by Solomon [9] and Haigh and Skelton [10] are shown in Figure 3. 

The scatter band of solomon's correlation is a factor of slightly more 
than 3, and the scatter band is reduced to a factor of 2.5 by the Haigh 
and Skelton correlation. 

Brown, et al. [11] calculated the plane stress crack tip field of 
Solomon's specimen with FEM. In their calculation, a constant stress 
boundary condition was used. The crack growth rate was correlated with 
the size of the severe strain zone, r g . The constant stress boundary 
condition will introduce a bending moment at the test section, and the 
bending moment could be more than that experienced by the specimen. 

Solomon's specimens were tested in the strain controlled fatigue test 
fixture. Both ends of the specimen were firmly attached to a rigid test 



5 


frame. Therefore, the test condition can best be simulated by a constant- 
displacement boundary condition, which was used in the present FEM calcu- 
lation. 

For a specimen under a general-yielding cyclic-load, the entire test 
section of the specimen experiences cyclic plastic deformation. The applied 
cyclic stress and cyclic strain as well as the crack tip cyclic stress and 
strain fields are related to each other through the cyclic stress-strain 
curve . 

The cyclic yield stress (Ao y /2 = 170 MPa) was obtained as the inter- 
section of two straight line segments in the log Ao/2 versus log Ae/2 plot, 
one in the clastic region and one in the plastic region (Figure 2b). The 
cyclic-s tress cyclic-plastic-strain relation derived from the 

data in Figure 2b is 

A 0.26 
At 

y (MPa) = 1360 (y 2 ) (9) 

where Ao/2 and Ae^/2 are cyclic stress amplitude and cyclic plastic strain 
amplitude, respectively. 

The complete cyclic stress-strain history experienced by the material 
in a crack tip region is rather complicated. As a crack tip advances 
ahead toward a point, the material at the point experiences increasing 
stress and strain ranges. As the cyclic loading goes on, the material at 
the point either cyclically hardens or softens. The detailed 
cyclic constitutive relation of the steel is unknown. Crack 
closure even complicates the matter further. In view of the 
complications of the crack tip field, it is necessary to make a 
few simplifying assumptions in order to calculate crack tip 
field for fatigue crack growth analysis ^ 
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During the loading half cycle, the material near" the crack tip in 
the plane stress region near the specimen surface will be stretched. 

Upon unloading, the crack front in the plane stress region will close. 

I : 

Crack closure will reduce the effective stress intensity range. It is 
well known that the compressive residual stress and crack closure due to 
a tensile overload will cause crack growth retardation. Normally fatigue 
crack growth is measured under a tension-tension load. If a tensile 
overload is followed by a compressive overload, the crack retardation is 
greatly reduced because the crack tip compressive residual stress field 
and crack closure are washed out by the compressive overload [13]. For 
the completely reversed loading of low cycle fatigue, the large compres- 
sive strain will "flatten" the crack surface and will even up the residual 
stresses throughout the specimen. Therefore it is reasonable to assume 
that at the beginning of the tensile half cycle, the specimen is stress 
free, the crack tip will open as the applied stress becomes tensile as 
suggested by Haigh and Skelton [10], and the fatigue crack growth mechanism 
will become operational. 

We also assume that the material at a crack tip experiences "stabili- 
zed" hysteresis loop. This assumption is reasonable if crack growth rate 
is slow enough and the material has experienced enough number of high 
strain cycles. In the small scale yielding case, the material will have 
to experience several thousand fatigue cycles in the monotonic plastic 
zone and, in addition, it will have to experience several thousand cycles 
of cyclic plastic deformation before, it reaches the crack tip. Therefore, 
the hysteresis loop of the material close to the crack tip is fully 
stabilized. In the case of low cycle fatigue, the number of fatigue cycles 
is much less, but the cyclic strain range is much higher. Because of the 
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high strain range, it is also reasonable to expect stabilized hysteresis 
loop. If hysteresis loops are stabilized, the cyclic stress-strain rela- 
tion can be used for the finite element calculation. One will be able 
to calculate the cyclic stress and cyclic strain fields by making a 
static calculation using the cyclic stress-strain relation, if the hys- 
teresis loop is stabilized everywhere, and if the residual stress is 
negligible at the beginning of the tensile half cycle. 

For Solomon's specimens, the size of the plastic zone is much larger 
than the plate thickness, therefore, plane stress finite element should 
be used. The mesh is shown in Figure 4a. The lateral curvature of the 
specimen was simulated by layers of different thicknesses, as illustrated 
in Figure 4b. The detailed mesh near the crack tip is in Figure 4c. The 
solid lines in 4c delineate the meshes, and the dashed curves are the 
paths for J-integral evaluation. All the elements in the layer closest 
to the crack line are of the same thickness, and the values of J-integral 
are evaluated along the paths within this layer. 

The ABAQUS FEM program was used. The plane eight noded quadrilateral 
quadratic isoparametric element was used. The crack tip triangular elements 
were formed by collapsing one side of the quadrilateral element. 

According to Hutchinson [14], and Rice and Rosengren [15], the crack 
tip stress and strain fields for non-linear elastic solids obeying the 
power-law stress-strain relation, (o/o q ) = a(e/e o ) n , are 



o 

o 



n 

1+n 

‘o'. . (6,n) 
ij 



e 

o 



(10a) 


(10b) 
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According to Equation 10, J is capable of characterizing the entire 

crack tip stress and strain fields. However, the recent' plane-strain 

finite element calculations (16, 17, 18) indicate that crack tip field ■ 

characteristics are strongly dependent on specimen geometry as well, as 

load level. In other words, at the same J-value, the crack tip field may 

vary widely from one specimen geometry to another. For the same specimen 

shape, at the same J-value, the crack tip field in a small specimen in 

deep general-yielding may differ considerably from that of a large specimen 

at a limited amount of plastic deformation. According to Equation 10, if 

J is capable of characterizing crack tip field, the plot of (c^/o^) or 

(c../e ) versus r/(J/a ) should fall on the same curve regardless of the 
ij ° o' 

spceimen geometry and load level. . 

Crack tip stresses and strains for three crack lengths, a = 1.27, 

2.54 and 5.08 mm, were calculated. The data on effective stress, ( 0 ^/ 0 ^) 
and effective plastic strain (e^/Cy) for a = 1.27 mm at various load levels 
are plotted against r/(j/Oy) in Figure 5a. All the data in the crack tip 
region fall on the same curve. The slopes of the lines close to the crack 
tip are 0.85 and 0.23 which are slightly higher than the values given by 
Hutchinson [13] and Rice and Rosengren [14]; 1/(1 + n') = 0.8 and 
n'/(l + n') =0.2 respectively, n' is the cyclic strain hardening exponent 
Similar plots for o and e are shown in Figure 5b. Figure 6 shows the 

K yy yy 

normalized strain distribution, ( r. /ty) versus r/(J/oy), near the crack 
tip for three different crack lengths. In the crack tip region, the data 
for all three crack lengths fall on the same straight line as suggested by 
Equation (10). The data in figures 5a, 5b and 6 indicate that J-integral 
is capable of characterizing the crack tip field. In other words, at 
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the same J-value, the crack tip stress and strain fields are the same 
regardless of the length of the crack. Since crack growth rate is con- 
trolled by crack tip field, J will be able to correlate well with Solomon's 
crack growth data. 

Solomon's crack growth data are shown in Figure 3. At each of the 

data points, one can obtain the values of da/dN, At^ and crack length, a. 

AK = {(Ac /2) + At } /na . With Ae and At known, the value of the 

cep e p 

crack length can be found. With both crack length and the plastic strain 
at x = 7.62 mm given, the value of J can be calculated. The results of 
the calculations for three different crack lengths, a = 1.27, 2.54, and 
5.08 mm, are shown in Figure 7. The values of J are plotted against 

both c and e p at x = 7.62 mm, the location of the extensometer . 

yy yy 

The crack growth data were measured by Solomon at six different 
plastic strain levels, Ae p = 0.001, 0.002, 0.005, 0.01, 0.02 and 0.05. 

The J-values at these six strain levels were calculated for the three 
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In this part of the loading cycle, the relation between stress and 
strain follows the curve: ABC. This: segment of the stress-strain curve 
consists of the elastic and- plastic parts. The plastic part is Ae^ and 
the plastic part is h.z^/2. The value of AJ for Solomon's specimen was 
evaluated at the plastic strain range Ae^ tested, and it corresponds to 
the AJ along the loading curve ABC in. Figure 9. 

The stress-strain relation along ABC differs slightly from the 
cyclic stress-strain curve. However, the cyclic stress-strain curve was 
used to calculate the crack tip field for Solomon's test. The values 
of J are strongly dependent on strain and much less on stress, therefore 
this approximation will introduce a rather small error. Furthermore, 
the stress-strain relation along ABC depends on the applied strain range. 

It varied from one specimen to another and it is unknown. 

The correlation of da/dN with AJ for Solomon’s test is shown in 
Figure 8. The correlation is better than with either A(PK) or AK^ 
as shown in Figure 3. The data in Figure 8 give the empirical relation 

(mm/ cycle) =0.7 x AJ^’^ (11) 

dN 

for Solomon's 1018 steel. AJ is in MPa-m. Equation (11) is also plotted 
in Figure 1 as the dashed line. The data for 1018 steel are within the 

i 

data band of A533B steel. It is well known that fatigue crack growth 

data of low and medium carbon steels are all close together. 

However, in order to use the correlation for engineering applica- 
tions, it is desirable to compare the measured and calculated AJ. 

CONCLUSIONS 

In general-yielding cyclic-loading, fatigue crack growth rate cor- 
relates very well .with the calculated AJ values for specimens tested at 
six different cyclic strain ranges. In order to predict fatigue 
crack, growth rates in engineering structures, it is necessary 
to establish the correlation between da/dN and calculated AJ. 
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Figure 3. 
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Figure 5a. Normalized effective plastic strain (e^) and effective 
stress (o £ ) distributions ahead of the crack tip for 
Solomon's specimen. 
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Figure 7. J-integral dependence on applied strains at x = 7.62 mm, 

the location of the monitoring extensometer . 









Figure 8. 



Cyclic stress-strain loop and the deformation work 
density, W^. 
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